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IMPORTANT ANNOUNCEMENT 


The increase in the output of papers in the field of mathematical 
biophysics makes it difficult to insure prompt publication without an 
increase in the size of the journal. Therefore, the Bulletin of Mathe- 
matical Biophysics inaugurates the following service: 

Upon acceptance of a paper, the Editor, if necessary, will ask the 
author to shorten the paper to an extent dictated by the requirements 
of a reasonably prompt publication. The shortening should in no case 
reduce the paper to a mere abstract. Such a shortened paper will be 
published within six months or less. 

The unabbreviated original manuscript will be kept on file at 
the editorial office. Any person desiring to avail himself of the com- 
plete manuscript, may obtain promptly a microfilm copy of the latter, 
at the cost of 1¢ per page plus postage, by applying to the Editorial 
Office, 5822 Drexel Avenue, Chicago, Illinois. 

All papers in the Bulletin which have been thus shortened, will 
be marked at the end by the symbol MF, followed by a figure, indicat- 
ing the number of doublespaced typewritten pages of the unabbrevi- 
ated manuscript. 
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SOME THEORETICAL CONSIDERATIONS CONCERNING THE 
INTERCHANGE OF METABOLITES BETWEEN 
CAPILLARIES AND TISSUE 


INGRAM BLOCH 
THE UNIVERSITY OF CHICAGO 


The author’s previous treatment of diffusion of oxygen from an 
idealized capillary is made more exact; the treatment here is valid for 
any substance which diffuses between blood and tissue, though the physi- 
cal situation is still somewhat idealized. A general equation is derived 
but not solved. 

Then follows an approximate treatment of the case in which all 
capillaries within a certain sphere have ceased to flow; conditions are 
discussed under which cells at the center of this sphere will die. 


In an earlier paper (Bloch, 1941), the author presented an ap- 
proximate treatment of the problem of diffusion of oxygen from a 
capillary. The present paper contains a formulation of the exact prob- 
lem of diffusion of any substance from or into a capillary. 

We shall again consider a circular cylindrical capillary of length 
lL, whose axis is to be the ¢-axis, so that ¢ = 0 at the arteriole and 
€¢=1, at the venule. r(P) will be the distance of point P in the tissue 
from the ¢-axis; 7) will be the radius of the capillary, and v the ve- 
locity of the blood-flow. D and h will denote, respectively, the tissue 
diffusion coefficient and the capillary permeability of the substance 
under consideration, while Q(P) will be the rate of production of the 
substance in gm/cc seé at point P in the tissue. 7, and h are assumed 
to be constant. | 

I. The case of oxygen: For consideration of oxygen, we shall 
need to define some more quantities: 

c,(C) =concentration, in gm cm, of oxygen dissolved in the 

blood. 

C,=¢,(0), C2 = ¢(l), c'(¢) = concentration, in gm cm~, of 
oxygen in the tissue immediately outside the capillary 
wall. . 

x.(¢) = concentration, in gm per cm’ of blood, of removable oxy- 

gen in oxyhemoglobin. 

t1 = x(0), Xe = Xo(1), To(C) = 2(F) + & (fC); T,=%+¢,, 
T, = 2, + ¢,¢c(P) = concentration, in gm cm, of oxy- 
gen at point P in the tissue. 


1 


2 INTERCHANGE OF METABOLITES BETWEEN CAPILLARIES AND TISSUE 


2p ef ee (1) 
Ci mae ne A 


From material balance (Bloch, 1941), we have the following ex- 
pression: 


—vry = = 2h[e.(C) — ¢'(6)] = 8aDo. (2) 
Also we know that the diffusion equation may be used in the tissue: 
Q 
2 a 
Ve D: 


The diffusion equation may be otherwise stated; the capillaries in the 

volume V’ of tissue will make a contribution to ¢(P) which will be 

given by the sum of several surface integrals, each of which is taken 

over the surface of one capillary; added to this sum will be a volume 

integral, which represents the effect of tissue metabolism within vol- 
o (Ci) ds; 


ume V’: 
Bis ore 2, ecun yee lal [aren 


Here S; is the surface of the i-th capillary: p(pi,P) is the distance 
from P to p;, the general point on S;, and p'(P’,P) is the distance 
from P to P’, the general point in V’. Equation (3) is a special case 
of the general solution of Poisson’s equation (MacMillan, 1936). 

Bloch (1941) gave a brief discussion of the effect of other capil- 
laries on concentration of a metabolite at a point in the region of 
supply of some particular capillary. The only conclusion reached was 
that, in general, this effect could not be neglected. It will, however, 
be evident that distant capillaries will have less effect at points im- 
mediately outside the surface of some particular capillary than at 
points nearer the boundary of the region of supply. Hence, we shall 
set up equation (3) for c’(z), and shall neglect all terms but the first 
in the summation. The z-axis is to be taken coincident with the ¢-axis, 
but whereas ¢ is defined only for points within the capillary and on its 
surface, z is defined only for points in the tissue. Thus c’(¢) and ¢’ (z) 
are logically distinct, but physically identical. 

Equation (2), combined with equation (3) set up as proposed 
above, will still require one more relation to make the problem deter- 
minate. This relation may be most conveniently obtained by consid- 
eration of the equilibrium between dissolved oxygen and oxyhemo- 
globin, which consideration should lead us to an expression for x, (¢). 


INGRAM BLOCH 3 


%o(¢.) having been found, equation (3) for c’(z) can be solved, 
and the quantities c’(£), co(¢), 2(C), T,(¢), and o(C) will all be de- 
termined. Then this value of o(¢) can be substituted in each term of 
equation (3), which substitution, combined with a reasonable assump- 
tion about Q(P), will reduce the problem of determining c(P) in the 
field of n equivalent capillaries to that of evaluating the integrals. If, 
on the other hand, we prefer not to make any assumptions about 
Q(P), we can use H. D. Landahl’s expression for Q(c) (Rashevsky, 
1940) in the volume integral, in which case c(P) will be the solution 
of an integral equation. Thus, by the use only of the approximation 
involved in dropping all but one of the terms in the summation in 
equation (3) when solving for c’(z), we shall be able, in principle, to 
obtain an expression for c(P). 

Bloch (1941) gives a discussion of the equilibrium between oxy- 
gen and oxyhemoglobin, in which it is pointed out that x, will depend 
not only on ¢,, but also on concentration of CO, and on the time lag 
in dissociation of oxyhemoglobin. It has seemed expedient to neglect 
the time lag, and to determine the relationship between x, and c, and 
CO, concentration by reference to data in L. J. Henderson, (1928, p. 
130). The curves here show y , percentage of hemoglobin oxidized, as 
a function of w, pressure of dissolved oxygen in millimeters of mer- 
cury. The parameter of the family of curves is w, pressure of CO; 
since its effect is slight (if w and u are expressed in mm Hg (du/dw), 
= 1/4 per cent mm'(dy/dw), = 2/3 per cent mm), we have set 
y(u,w) = y(u, average w). y(u) may be fairly well represented in 
the significant range of u( from 20 to 80 mm) by a parabola: 


U——O5U I Mlb 1 ly 


where o. = — 2.18 X 107, o, = 3.24, and a = — 27.3. 
To get x.(¢,), one must get x,(y) and u(c,). Both these relation- 
ships are simple proportionalities. 


——! Us 
~ 100 


Xo Y= My (4) 
where x, is the concentration of oxygen in HbO.(gm cm* of blood) 
when Hb is 100% oxidized, and M = x,/100.* Also, u is proportional 
to Co: 


Uu=Ne. 
11It may be of interest here to note that anemia will be characterized by a 


low value of the quantity «,. Perhaps in the future it will be possible to inves- 
tigate anemia from this point of view. 
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Hence we have 
= My = M (anu? +0, + 0%) = M (a,N2%C,? + oNC. + Qo). (5) 


Since 75 = % + Co, 
lie — Co aia M (a,N?C,? ta a,NC, ola Oo) = Y20o" =e 10 == Yo (6) 


the y’s-being defined in terms of the a’s and N by equation (6). 
Now we are able to combine equations (2), (3) and (6) into a 


single equation in one unknown. 
First, since we shall be manipulating only the surface integral in 
equation (3), we shall substitute A for the volume integral. We have 


BORA car 


¢' (2) = [0 | eee pee +A (7) 
1/272 (licos OV (2 — 6)" 


where @ is the angle between the radii vectores to the points p and P. 
Only one angle appears in the integrand because o is assumed to be 
a function of ¢ only. 

From equation (2) we learn that 


or 


_ ro AT 
o(f) = BaD de (8) 
and that 
st al vr, AT, 
e'(f) =eo(¢) + 2h di” (9) 
Equation (6) gives us another expression for dT,/d¢: 
dT yas de 
dt = (2y26o + v1) rm (10) 


If we substitute in (8) and (9) the value of dT,/dt given in (10), 
and substitute the resulting expressions for c¢’(¢) and a(C) in (7), 
we get 


Co (Zz) Syed rae 5, [2260 (2)ee- n 


(gg fe ones de, 
7 co f SnD [2y2e0(f) + v1] “ae (11) 
S de rs Q(P’)dV’ 

\V/2r,?(1 — cos O) + (2 — CF)? pul Do (PEPYS 
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Equation (11) will, in principle, give us ¢(¢), from which we 
can readily obtain ¢’(¢), «(¢), T.(¢), and a (¢). If we put the value 
of o(¢) thus derived into equation (3), we will have a deterininate 
expression for c(P). Therefcre the present desideratum is the solu- 
tion of (11). Unfortunately, however, equation (11) is of a form 
which is, so far as is known, not soluble. Hence we shall have to con- 
fine ourselves to an incomplete discussion of the solution. 

First, it is evident that, since O S z S 1, there will, for every 
choice of P(z), be one point in the range of integration in (11) at 
which the integrand will become infinite. That this singularity will 
not make the integral become infinite can, however, be shown by the 
following argument: Consider a small circle, of radius 6, fitted to 
the surface of the capillary with its center at the point P(z). If 6 
is sufficiently small; we may approximate the integration over the in- 
terior, 7, of this circle by integration of a constant oc, the average 
value of o in the neighborhood of P, over a circular region of a plane; 
as 6 approaches zero, this approximation becomes exact; for any 6 > 0, 
the singularity in the integrand occurs within the circle, so that there 
will be no singular points encountered in integrating over such of the 
capillary surface S as lies outside 7. Now, if « denotes the distance 
from the center of the circle to the general point in 7, our integral 


becomes 
5 6 
O° Qe — _ 
[ i= i 2aodse = 2100. 
0 é 0 


This quantity obviously approaches zero as 6 approaches zero, so that, 
as 6 approaches zero, the integration over S — 7» approaches a finite 
limit; it is evident that the singularity at P does not make the integral 
diverge. 

II. The case of glucose: The problem of the diffusion of glucose 
is simpler than that of oxygen because glucose experiences no buffer- 
ing so, 


Co(f) =To(C), (tl) =0. 
For this reason, the equation corresponding to equation (11) is 


Be U —Uro 
ee oe ‘of 2 +A. (12) 
¢o(2) 2h dz : ‘ 8nD /2r,2 (1—cos O) + (z—C)? 


The same argument as that given for the case of oxygen applies here 
and shows that the singularity in the integrand does not make the in- 
tegral diverge. 


6 INTERCHANGE OF METABOLITES BETWEEN CAPILLARIES AND TISSUE 


III. Other cases: It will be seen that other metabolites will be- 
have either like oxygen or like glucose; we shall denote the types of ' 
substances exhibiting these two types of behavior as types I and II, 
respectively, type II being, as we have seen, a special case of type I. 

A. Non-buffered substances: All non-buffered substances behave 
like glucose, and hence belong to type II. Therefore equation (12) 
applies to their concentrations in the blood, and its solution will en- 
able us to compute all other unknown quantities. 

B. Buffered substances: In the general case of buffering, the 
substance will exist in several forms, in some of which, the ‘mobile’ 
forms, it can pass through the capillary wall and in the rest of which, 
the “immobile” forms, it cannot. Let us denote by c,.(¢) the total con- 
centration, in gm per cc of blood, of the substance in mobile form, by 
Xo(C) its total concentration in immobile form, and, as for oxygen, by 
T,(C¢) the sum ¢c)(C) + x4 (C); let c(P), x(P), and T(P) be the corre- 
sponding concentrations at point P in the tissue. Let us assume we 
know « as a function of e(2, as a function of ¢,) when the forms are 
at equilibrium. 

We are treating a steady state, so we assume the existence of 
equilibrium between mobile and immobile forms in the tissue, and, in 
accordance with the treatment of oxygen, neglect any time lag there 
may be in the transition between the two forms in the blood. 

It will be seen that buffering in the tissue will manifest itself by 
influencing the production integral in equation (3). But, since we are 
dealing with the case in which exactly as much of the substance goes 
from the mobile form to the immobile form as goes the other way, 
at any point in the tissue, it is evident that buffering in the tissue has 
no effect on the value of Q, and hence introduces no change in our 
problem. 

Buffering in the blood will manifest itself exactly as it did in our 
treatment of oxygen. o will be defined in the same way, and equations 
(2), (3), and (7) will be unchanged. The steps between equation 
(7) and equation (11) will be essentially the same as before, except 
that, in general, we shall not have x, as a quadratic function of ¢, . 

dT, dx deo : 


d 
dt ae [%(¢o) + Co] = - ai ee 


dx E 
a + 1] will be known; let us denote it by $(c). Hence, from equa- 
(0) 


tions (8) and (9) 
eee Vo dco 
RAND sheer sh AG Re (13) 
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e'(C) =¢,(¢) 1 oa (4) (14) 


The generalized version of equation (11) will thus be 


Co(z) + = = (9) 


3 d. 
Ur ,2 (Ca) —— a 
do | — — eo 
E ‘ 82D \/2r.?7(1 — cos ©) + (z—£)? 


It will be evident that equations (12) and (14), for general sub- 
stances of types II and I respectively, will hold whether the substances 
are being produced or consumed in the tissue, since the sign of Q does 
not affect the validity of equations (2), (3), or (7), from which the 
later equations are derived. Actually, also, equation (12) is a special 
case of equation (15), in which ¢(c,) =1. 

Thus the solution of equation (15) will determine the behavior of 
any metabolite in its passage between blood and tissue. 

IV. Since the solution of equation (15) presents seemingly in- 
surmountable difficulties, we shall derive an expression which should 
enable us to find, approximately, the effect of distant capillaries. 

First, because we are interested in the concentration at a point 
far from a capillary, we shall treat the capillary as a line segment 
of length /, omitting the integration with respect to © in equation 
(3). In accordance with this change, we must replace o, productivity 
per unit area, by 1, productivity per unit length. It will be seen that 


A=2arQo . (16) 


(15) 


Hence, if we denote by o;(P) the contribution to c(P) made by 
the i-th capillary, and by A(P) the effect at P of production in the 
tissue, as in equations (7), (12), (15), we have 


c(P) — 2 a(P) + A(P)- (17) 


The quantity c;(P) will be given by the following equation: 


Li 

A(C1) dei oo (G5) de 
HY gd a ———- = 2a” ——$_—______—_—-.. 18 
ase J p(6i,P) anf Vr? + (2: — 6)? mad 


- Determination of «(¢) requires the solution of equation (15) or 
an equivalent equation, in which can not be used the line-segment ap- 
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proximation of this section, since the equation applies to a point ad- 
joining the surface of the capillary. For this reason, we shall inves- 
tigate the form of c(P) for two simple forms of o , assumed given. 

From equation (2) we note that, if «(¢) is a polynomial of degree 
K, T,(C) is a polynomial of degree K + 1. In particular, if o(¢) is 
a constant, o(£) =o, we have JT,(¢) a linear function of ¢, 


of vod Es 
TG are eee 
and eee (19) 
Vie eA) weeny 


Also, if « is linear, T, is quadratic: 
T.(C) = To aif bal oF bel? 


and a (20) 
o(f) =— BaD (ur + 2u2C) 


In the case of quadratic T,(¢) [equation (20)], we have not 
enough facts at our disposal to evaluate the coefficients u, and we , but 
if T,(¢) is linear [equation (19)], all coefficients are determined. It 
will be seen later that this linear form is the most useful. 

Us Lar ls 
82D l : 
Thus, from equation (18), we have, dropping the subscripts and per- 
forming the integration, 
2 a 2 = 2 = 
c(7, z) — re (Ps — Ta) log MARES Ded eric (21) 
ADI . Vee 2 


If T, is quadratic and o linear, as in equation (20), so that 


If «(¢) is constant, equation (19) tells us its value is 


Ur. 


o(L) utah a7) (ta a 2ucl) ? 


then integration gives 


Wee 1 HR eat eee eae 
(1) #) =) cab te a ee 
r+ (e—N?t+l— (22) 
ide ht) aloe ee 
WTR 48 


When one is dealing with an actual capillary, T, and o cannot 
be linear and constant, respectively, since expression (21) is sym- 
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metric about the plane z = I1/z, while linearity of 7, implies asym- 
metry in c, and in c’, hence in e(P), if permeability, h , is assumed 
constant. However, in treating the case of many capillaries, oriented at 
random, we shall assume that for each capillary with a given orienta- 
tion there is another at the same distance with opposite orientation, 
and shall use equation (21) to express the effect of one capillary. 
Distance in this case will be the distance from the central point of 


the capillary: 
S, l 
R=,/r ar tao hee 
z 


or, if 


2--=Z,R=VPF2. 


aQ 1s 


We shall set r= Z on the average. Now equation (21) becomes 


= ries sl a). 
e(R)( = c(7, Z) Pare 


Nee (23) 


Since we are interested only in the effects of distant capillaries, 
we shall assume / is small compared with R and shall hence subtract 
4/2 from the radicands in both numerator and denominator of equa- 
tion (23), thus making the radicands perfect squares. Now it is easy 
to integrate c(R) over the volume between R = R, and k = R,. This 
integration will give us an expression for contributions to concentra- 
tion at point P of all capillaries lying between two spherical surfaces 
of radii R, and R,, respectively, with their centers at P. Rh, > F,. 
What we are doing here is treating these capillaries as a producing/ 
consuming continuum rather than as discrete producing/consuming 
bodies; we shall obtain } c:(P) [see equation (17 )] for the case in 
which all capillaries within a sphere of radius R, have ceased to flow. 


if the tissue has radius FR, . 


10 INTERCHANGE OF METABOLITES BETWEEN CAPILLARIES AND TISSUE 


+. herd mnvr,?(T wo T>) Pp R. ai 41 
208) STs) eal eas cre 
(24) 
Padua a R, +41 
+ = log ate oe ae + 41(R.? — Ry2){ 
RY-— i a 
4 


Here n is the number of capillaries per cubic centimeter between R, 
and R,. A(P) will be given by the following integration: 


Re 


= Ole tani se Ore 
AP) = 75 | “a R= ap Re. (25) 


Thus, by combining equations (24) and (25) in accordance with 
equation (17), we obtain 


mnvr,? (Ts, — ES) R, + 4l R,+4l 
pp re eS 2s 3 
c(P) SDI E log eas R,? log Rees 
2 
oe ee 
Ih a ces 2 eG OY 
xe pee cal Ske | +2 ne 
ji bs 


4 


Equation (26), then, enabling us, as it does, to calculate concen- 
tration in a region whose capillaries are not flowing, should be sus- 
ceptible to experimental verification. 

An assumption which must be made concerning tissue is that, if 
equivalent capillaries are uniformly distributed, m per unit volume, 
throughout tissue, each capillary produces/consumes exactly one n-th 
as much of any metabolite as is consumed/produced in a unit volume 
of tissue. This must be true if the tissue is in a steady state. 

If, as in the case we are considering, the capillaries inside a 
sphere of radius R, are not flowing, then capillaries at distances from 
the center not much greater than R, will have to supply/consume 
some of the material that normally would have been supplied/con- 
sumed by the capillaries inside the sphere. For this reason they will 
have values of o greater in absolute magnitude than the average value 
for the whole tissue; hence equation (26) is not accurate, since o in 
fact should be a function of R and so should have been included in the 
integrand. This correction we shall neglect. 
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Whatever may be the correction to o for the closer capillaries, it 
will be seen that the more distant capillaries will be less affected by 
the absence of capillaries within the sphere of radius R,; as R ap- 
proaches infinity, the correction approaches zero, so that capillaries 
at infinite distance from P will each produce/consume exactly as much 
material as is consumed/produced in one n-th unit volume of tissue 
at infinite distance. Hence, 


so) dc(P) 
lim = 


R20 d 


(27) 
2 

Application of equation (27) to equation (26) will give us a re- 
lation between Q and (7, — T.) in terms of the other quantities in- 
volved: 


l 
R,+- 
RNvT 3? (T, — To) ey Q_ 
py EEO pt H#),9=o, (28) 
aes? 
whence 
mnvr,?(T, — T2) Ampere 
deat fA nee TD lea =0 
or 
nae 29 
nur, (T ss T's) a ae me ( ) 


Equation (29) is identical with the relation between Q and 
(T, — T.) obtained by a different method. 

If blood is flowing with velocity v in a capillary of radius 7 , then 
the volume of blood passing a given point per unit time is 27’v . Then, 
if the concentration of some substance is 7, at the entering end of the 
capillary and 7. at the other end, the blood entering per unit time con- 
tains a mass of the substance equal to T,17,?v, while the mass con- 
tained in the blood leaving per unit time is T.77,?v. Hence, in unit 
time, the blood loses (JT, — T2)x7?v, or gains (T; — T,)ar.?v, mass 
units. Now, if there are n capillaries per unit volume of tissue, and if 
we assume each capillary to produce/consume as much of the sub- 
stance as is consumed/produced by a volume 1/n of tissue, then 


(T, — T2) aru = — £ (30) 
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if Q is production per unit volume of tissue. Equation (30) is seen 
to be identical with equation (29). This fact does not constitute a 
verification of our assumption that o is constant, since the assump- 
tions leading to equations (29) and (30) are equivalent for any a(C), 
and thus must lead to the same equation. Hence the identity of (29) 
and (30) merely indicates that the error introduced in the approxi- 
mation used in integrating equation (23) vanishes as R approaches 
infinity. 

If we now substitute for Q in equation (26) its value as given by 
equation (29), we obtain 


abe mvr,?(T x T2) 5 Re si 41 ‘ Tie =F 41 
c(P) ST Ty Te R.' log Real R,} log ee 
Bae Ra tae ees ee oe (31) 
a 3 log Ba th 5 ies 5 IR, 


It is now possible to obtain an expression for concentration at 
any point P at distance R from point P, if R < R, and if we assume 
the field within the sphere of radius RF, about P is spherically sym- 
metrical. This comes from solution of the diffusion equation 


Q 
oe 
Wie 5 
whose solution for this case is 
ANT»? (T, — To) R.+ 4 
ee Ae ee eR ee 3 
c(R) ani E ] at 
Ri+#l PF. R2~-4P : (32) 
— R? log — ope eS Ba 3 Si a pain 
og 1 — 4 oF 8 loge 9 iP 41R,? 5 IR. =F 4IR | 


All the foregoing treatment has been carried through on the as- 
sumption that Q is constant. Whether or not this will in general be 
true is difficult to say; Landahl has, however (Rashevsky, 1940), ob- 
tained relations between concentration and consumption for glucose 
and oxygen, respectively. These relations show that, as concentration 
of either of these substances outside a cell increases, the cell’s con- 
sumption of this substance approaches an asymptotic value —Q*. It 
is possible to define in each case a value c* of external concentration 
such that Q may be considered equal to Q* if and only tives co 5it 
seems also reasonable to assume that metabolism of other substances 
depends upon the metabolism of oxygen and glucose in such a way 
that if oxygen and glucose concentrations, denoted by C and C’ re- 
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spectively, are greater than their critical values, C* and GC”, all produc- 
tions and consumptions are constant, provided the constants of the 
cells in the tissue do not change. Hence, our treatment is valid if 
and only if C 2 C* and C’ 2 C™. In Rashevsky, 1940, we find expres- 
sions for GC and C’: 


C=2% Ot GO (33) 
eee EQ’ 
C — Xx QO” a5 oO: =e Q’ S (34) 


Here Q, Q*; Q’, Q* are consumption and critical consumption of oxy- 
gen and glucose, respectively. €, &’, y, y’ are constants given in Ra- 
shevsky, 1940; my zy is Rashevsky’s ¢. These constants depend, in par- 
ticular, upon the nature of the cells in the tissue. 

Now we set 


Gray t+é 


35 

Ge — x a br. ( ) 

These are the values of C and C’ at which the straight lines tangent to 
the C(Q) and G’(Q’) curves, respectively, at the origin, intersect the 
lines Q = QO* and Q’ = Q”. They seem reasonable values for C* and C”. 
If we now substitute for T., T., D the values T,, T. 2, Tye 

T,', T.', D’ for the special cases of oxygen and glucose, respectively, 
and substitute for c(P) in equation (31) the values of C* and CG” from 


equation (35), we obtain 


MNUT 2? (T1 — Te) R, + 4l 
[591 Sah SS ma Be tO es 1 


reo wnt eameaad ba ee 3 
ae Tess oe ayy ats 8 ‘Cay area? = 41R{? BE 3 lR.? 
(36) 


yt F&= 


nNnvr2(T1' — T2') ee R,+ 41 
8D 


These equations give the conditions that concentrations of oxygen 
and glucose, respectively, are as low as they can be if metabolism is to 
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remain constant as concentration changes. If R,* and R,”* are the 
values of R, for which these two equations are true, then the smaller 
of these is the radius of the largest sphere of capillaries whose ceas- 
ing to flow will not decrease metabolic rate at any point of the tissue. 

If we divide each equation (36) by the (positive) coefficient in 
front of the brackets on the right side, we get two equations, in each 
of which the quantity on the right side is the same decreasing func- 
tion of R,. Hence, R,* < Rk,” if 


3D1(xy + €) BDI Fe) 
anvT.- (TT, — Ts) ANUL C24 a a) 


3% 
or if (37) 


Dat) ee aaa) 
TT labia Te aT 


Thus we see that if the inequality (37) is true, oxygen is the critical 
substance in maintenance of asymptotic metabolism; otherwise, glu- 
cose is the critical substance. 

Once we have discovered whether or not inequality (37) is true, 
we know which of the equations (36) to solve for R,. It is to be ex- 
pected that when R, exceeds this value, the metabolism of the cells in 
the center of the sphere of non-flowing capillaries will considerably 
slow down, and in some tissues, notably in brain and in striated 
muscle, the cells may begin to die. 

The author wishes to acknowledge the helpful advice given him 
by Dr. H. D. Landahl, of the University of Chicago. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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The frequency of periodic elongations of a cell has been calculated 
as q function of the cell volume for Rashevsky’s dynamical theory of 
a cell. 


In the interest of more adequate experimental verifications of the 
existing theory of cellular pulsations, it becomes desirable to make 
some sort of calculations of the frequencies, as given by the theory. 
While any exact or absolute value of the frequency cannot be made 
because of the present lack of experimental data, it is possible to de- 
termine how the frequency should vary with the various factors in- 
volved. The following calculation consists of a discussion of the ex- 
pected variation of frequency as a function of the total volume of the 
cell, since it seems likely that variations of the volume may be per- 
formed most simply in an experiment. 

One might expect these calculations to be applicable to certain 
experiments, one of which might be amoeboid movements. An amoeba 
is, of course, highly irregular, but by measuring its area in a photo- 
graph, and dividing by its perimeter, some quantity which might be 
considered a rough approximation of its elongation would be obtained. 
This function might fluctuate very sporadically with time, but some 
approximate value of the average frequency would be found. The 
amoeba’s volume could then be changed by micro-dissection and thus, 
the curve of the variations of frequency with volume measured. 

N. Rashevsky (1939) has developed a theory of cellular dynamics 
where he finds 


d — —_— =e 
7 (a ~ 2) = Arla G2) + B,(¢ — C2) 


(1) 
ie 


(¢ — C2) = A2(a@— a) + B,(c — Ce) 
dt 


where a is one half the length of an approximately cylindrical cell, c 
is a linear function of the average cell concentration, and 
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A= -1)—-; A=--=|=—— 
( 2V. 270o%* lg —2 Jp. 
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w= (Fe) } oe Lacs Daler 


with V the volume of the cell, and the remaining notation the same as 
that used by N. Rashevsky. It was shown by him that undamped pe- 
riodic elongations occur if 


A, + B.==0; A,B, — B,A.>0. (3) 

Damped periodic elongations will occur if the first of these rela- 
tions fails, provided that the frequency w , given by 

a? AN Esa ese B.A; ey 1(A, =F Bove (4) 


does not have imaginary values. The vanishing of this expression, 
then, gives limits to the values of the parameters for which the cell 
may oscillate. It is the purpose of this paper to find the frequency as 
a function of the volume V of the cell. 


ol 122MD; 
na 2 [nod a) vn ai 


Set 
en ee 
V ‘ ~ 2n\86nD.2D; :  QayD;" (5) 
Substitution of equations (2) and (5) into (4) gives 
2 US ey ae 
vo? = y| —— — —— 
@ VWeses nTaleee 7 So 7 ) | : (6) 


N. Rashevsky has shown that a necessary condition for undamped 
periodicity is that 


Osis eV (7) 


and since the cell must produce (i.e. g > 0) to be able to oscillate, nec- 
essarily 


y=0 (8) 
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irrespective of the values of u or v, both of which are positive. 
Also, for real w , 


2 A* 1/3 
us E ae! J=°. (9) 


yt+3  4y eae U 


Equation (6) may be considered as defining vw? as a function of two 
variables, u and y. The equality sign will give the boundary of the 
region of the uwy-plane where equation (6) has positive values. Sim- 
plifications in this operation gives a quadratic in wu, 


Pee vgs Yo oN tickets Lae 
icc) Ay?/8 (ars apa te eal 
which has the solutions 
+ 
fe ee Oe (11) 
(Vy + V2)? 
epee wee (12) 
(Vy — V2)? 


It is clear that uw’, > u, for any y > 0 and consequently frequencies 
are obtainable over a region of the uy plane. It can be shown analyti- 
cally that wu, is a monotonic increasing function of y which starts at 
the origin with a vertical tangent and for large values approaches the 
curve u = y? from below. The wz’, curve rises from the origin (also 
with a vertical tangent) monotonically to infinity at y = 2, returns 
from infinity for y > 2 to some minimum value (wv, = 6.2) for 
y = 44.2, and then rises again approaching the curve u = y? from 
above for large values of y. For greater clarity, curves are presented 
in place of the detailed analysis. 

A simple numerical substitution will show that equation (6) has 
possible solutions for any given value of y provided 


UU es (18) 


This in turn implies that there exist solutions for 0 < wu < o with 
proper selection of y and conversely. As an example of this, the value 
y = 2 may be selected. Substitution in equation (6) gives 


Sea (14) 
5 W4Au 


which is a monotonically increasing function of wu for u > 0. For real 
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Ae a 
solutions of w, u must be greater than 5 y/4 but for large values, vw? 


4 : 
approaches the value 5 asymptotically. 


If one assumes a certain constant value for y while u varies, then 
vo? will be a function of u. The salient features of this relationship 
can be understood by considering vertical sections of the surface rep- 
resenting equation (6). The maxima of this family of curves of inter- 
section will be obtained by equating the partial derivative to ‘zero. 
Hence from equation (6) 


¥,) Vv 2 2/3 2 coast 1/3 1/3 2 = 
ae ee ere (15) 
ou A af AS U U Ves: 
Similarly by considering u as constant, 
d(var) ae eee ey “ay Eas y, 
oy yt3 (y+)? aula a 


(16) 
23 | oer oe 5 1 
ee ey | Sigs See eee | ee 
(5 eee ns U Tae ae te 


The first of these gives the locus of the extrema of the u-sections 
and has solutions 


y+3 : y+3 


ot te peas (17) 


when 0 = y < 2 and y > 2, respectively, u being always positive. wu, 
is a monotone increasing function, and uw’, decreases from infinity for 
y = 2 to a minimum and finally approaches the function y’’* for large 
values of y. 

The locus of the extrema for the y-sections, hich is of primary 
interest in this paper, follows from equation (16). Simplifying, 


(y = 2) (6 = 16y — 97) (a \ 2 y+6yt+6/ u 
So ARTE oe SRE of oO eee ie le ee es) 

(y + 3) y (y+ 3)? \y ; 
The curves represented by this expression are very complicated to 
discuss analytically and the results are shown graphically. There are 
in all six distinct parts to the set, two of which represent minima. 
The solutions may be found to be 
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= 
Us == Yy 2 (y + x 
3 oY (y + 8) (19) 


y+ 6y+ 6 =v (y+ 6y +6)? —8 (y—2) (y +38) (y+ 16y—6) 
(y — 2) (y? + 16y — 6) ae 


3 
u', == y8(y +3) X 


y+ 6y +6+4/ (y! + 6y + 6)? == (y= 2) (y +3) (y? + 16y — 6) 
(y — 2) (y* + 16y — 6) 


The upper expression is positive for all values of y and monotone in- 
creasing. For large y , it approaches the value y’’. The second equa- 
tion gives positive wv’; for 0 = y = 0.867 and 2 < y < o but for large 
y approaches the value 2y’. The first set is outside of the region of 
definition of equation (6) and denotes a minimum. This is true also 
for the second set for large values of y, but. between the value 2 and 
the minimum point, it represents a true maximum. These results are 
all shown in Figure 1 which is self-explanatory. 


0.1 ! 10 oo 1000 
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In order to study the volumetric dependence of the frequency of 
an oscillating cell, the analysis given above can be conveniently used 
to eliminate some calculations. In fact, for a given value of u (to be 
considered as a parameter) the expected variation of vw? as a func- 
tion of y may be sketched approximately, knowing the location of the 
zero points and the maximum. A single calculation of the actual value 
at the maximum will then give a clear idea of the expected relations 
which may be compared with experiments. These are shown in Fig- 
ure 2. 


FIGURE 2 


Preliminary computations of the parameter wu indicate that it has 
a very high value ~ 500, but the experimental data are very meagre 
and inaccurate, and much smaller values cannot be excluded from con- 
siderations. It is apparent that for u > 6.3 approximately, there are 
two separate regions for y in which oscillations are possible. From 
equation (5) large values of y imply small values of V, the cell vol- 
ume, provided V, remains constant. This would indicate that an oscil- 
lating cell would, as its volume is decreased, reach a quiescent state 
and then begin to oscillate again. It is entirely possible that this is a 
mathematical property of the equations without biological signifi- 
cance, but there seems to be nothing in the theory to preclude its pos- 
sibility, except the small size. 
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If we insert u = y¥?, the equation of the maximum for large y, 
into equation (6), 


cl ee Be 
vo =rheesc, , (yo) (21) 


which has the limiting value vw? = 


Along the curves u = y’? and y = 2, for very large u, both equa- 
tions (17) and (19) approach each other. Consequently in the region 
of definition, both partial derivatives vanish and we would expect ab- 
solute maxima for vw?. These values are given in equations (14) and 
(21) and it is consequently apparent that if v is a constant, there is a 
maximum attainable frequency of oscillation. 

The writer wishes to thank Dr. A. S. Householder for his careful 
guidance in the calculations and for checking the results, and Dr. N. 
Rashevsky for suggesting the problem as well as giving valuable bio- 
logical interpretations to the possible experimental implications. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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The development of a general theory of neuron-networks is here ex- 
tended to cases of non-steady state activity. Conditions for stability and 
neutrality of an equilibrium point are set up, and the possible functions 
representing the variation of excitation over time are enumerated. The 
inverse network problem is considered—which is, given a preassigned 
pattern of activity over time, to construct when possible a neuron-net- 
work having this pattern. Finally, a canonical form for neuron net- 
works is derived, in the sense of a network of a certain special topo- 
logical structure which is equivalent in activity characteristics to any 
given network. 


In essaying a treatment of the dynamical case, we shall find it 
coactive to take into account the finitude of conduction time. We may 
characterize each fiber of N by a total conduction time 6, defined as 
the sum of its proper conduction time and the average synaptic delay 
at the prevening synapse; it is clear that without loss of generality 
these quantities may be supposed equal. For at least one cannot redar- 
gue that their ratios are rational; and this being the case, they have, 
when expressed in terms of the smallest 6, a least common denomina- 
tor v. Replace, then, every fiber of total conduction time @ by a chain 
of v6 fictitious fibers and synapses, each of which has a total conduc- 
tion time 1/v and suitable thresholds, etc., so as to be otherwise the 
same as the replaced fiber. If this has been done, we shall measure 
time in units of 1/v in length, so that all total conduction times be- 
come unity. 

After these preparations we may write a set of equations for the 
y; consimilar to the set (2). The principal difference will be this: that 
if we view the quantities y; as functions of time, the excitation re- 
ceived by a given synapse s; from a chain c, in complete activity will 
no longer be determined by the contemporary value of the excitation 
at the afferent end of the chain, but rather by its value there for a 
time-point precedent by an amount equal to the total conduction time 
of the chain. In our units this latter is »;. We may therefore, write 
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y= ait DS os mi + 3S pix aj Bj As Ye (t — 5) 
; (1) 
Ag LA; iaeeca) +¥;(1 — 8:1} : 


If we use Boole’s operator E , defined by Ef(x) = f(x + 1) and 
having the obvious property E”f(x) = f(x + m), and change the 
origin for ¢t by setting t= t’ + «, where « is the largest of the »; , we 
may write this in the form 


Et y,= oi + gail yj +S pix Aj a; 8; BP? yn 
j i 


re (2) 
+ Ajlay(1— a) + (1 BI}, 
wherein we have set p; = ¢ — »; and dropped the prime on t. 
If we define the matrix 
H(E) = ||> $31 pix Aj 0; Bj BP4||, 
Jj 
this becomes 
[he — HE) ly ks (3) 


where RF is defined as before. It will be noted that when we give to E 
the value unity, the matrix H(F) reduces to M: this is consonant 
with the definition of the steady state as a condition where there is no 
change with time, which is to say, 


yi(t +1) = Hy (t) =yi(t) =1yi (2). 


If the matrix E* I — H be regarded as a polynomial matrix in the 
indeterminates EF , a; , 8; , we may, by the use of Smiths’ process, find 
unimodular matrices P and Q such that 


P[E*¢I—H|Q=T, 


say, is a diagonal matrix whose non-vanishing elements are the in- 
variant factors of E* I — H; let us denote these factors by 


T(E >a, B] ’ 


where i varies from unity to the rank of E* I — H and 7; divides T 
as is well-known. If we make the substitutions 


z2=Q'y,S=PR, 


equation (3) may be written 


i+1 9 


Tz=S, (4) 
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which, in scalar form, is 
Ti[E ,a, fp] 2=S; Bean lati ad beget (5) 


The system (5) is a set of independent difference equations in 
the z; whose coefficients are constants or multilinear‘in the multi- 
pliers. For any region I’, , these quantities will have determinate val- 
UES, Ap1, %p2, and the equations (5) may consequently be. solved by 
the standard methods, the solution holding throughout iis: It may be 
obtained explicitly as follows. 

We first derive the solution of the corresponding homogeneous 
equation, found upon replacing the right-hand side of equation (5) 
by zero. If we equate the function T;[E', a, 6] to be zero and solve 
for EF as a numerical magnitude, we shall obtain a system of roots, 
hia, di2 » +++, Ais , respectively of multiplicity 7, , no,---, ys, Say. Then 
the esaton in question is given by eae 


8 Ns iG 
SD D Gin 2's;, a: (6) 


j=1 k=0 


and the quantities a;;, are either constants or, more generally, arbi- 
trary periodic functions of period unity. 

To derive a particular solution for equation (5) we must distin- 
guish two cases. In the first, none of the roots 4;; is unity: we find 
easily by substitution that in this case a constant value.for z; satis- 
fies equation (5) ; this constant value is S;/T[1; a, 6]. In the contrary 
case there is somewhat greater difficulty. If unity be a root of 
T,[E ,a, 6], of multiplicity 7; , say, we may then put: 


TLE ,0, 6) 2= (E-1)" B(B)z.=A" B(B)aS Si, (7) 
where A is the difference-operator; if we make the substitution 
Vi= AZ, 
this becomes | 
P(E) u=Si, . (8) 


in which substitution of a constant value for v; is iden PeRIRES and, 
as before, yields v; = Si/W(1) as a particular solution... To find z;, 


we now have the equation 
A 2; =S;/¥ (1) 7 
which can be solved immediately as 


Si r(é+1) 


2 (1) Pe+1—n) 2G: +1)’ o 


y 
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This may be added to the solution of the homogeneous equation to 
yield the general solution of equation (5), which is accordingly 


Bs Sin out, VEE Ds Se § Ne-1 i i. os 
oP LY P(é+1—7)£5+1) FE Thin hd ije ( ) 
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“a 


The par ameters ;;x are to be determined from the particular circum- 
stances attending the entry of the network NX into the region in ques- 
tion. 

It will be instructive to compare the asymptotic behavior of thé 
solution (6), with the results of the purely static analysis made above. 
We note first i in this connection that the presence of unity as a simple 
or multiple root of some 7; is equivalent to the vanishing of the de- 
terminant | — M|:this follows at once from the fact that E*°I — H(E) 
becomes I — M when we set FE = 1, that T;[1,a, 6] are consequently 
the invariant factors of J — M, and that the number of vanishing in- 
variant factors of a matrix is equal to its nullity, which by hypothesis 
is at least unity in the present case. Considering first the case where 
(I — M| #0, we see that equation (10) assumes the form 


Si 
oe one 3 da Ave. (11) 


The 1;;, none being unity, may be divided here into three groups. 
First, there are those which exceed unity in absolute value: these con- 
stitute say the set ©, . Second, those such that |/;;| < 1; these may be 
collected into @,. Finally, those which are —1: these comprise 0;. 
Now if 0, and 9; are both null, the transient term on the left of equa- 
tion (11) tends to zero with ¢ independently of the initial values, so 
that the set: of y; correlated to 


2 —S,/Ti[1,a, p] 5 


if still within the region J’, in question, are the asymptotic values ap- 
proached by the network whatever its initial circumstances. We may 
therefore calf this equilibrium point a stable one. If this set of values 
is not within the region—which means that the multiplier-distribution 
corresponding to it does not satisfy the inequalities (8)—the network 
will always leave Ip) , however it enters. Now, if ; is null, but 0, is 
not, then in general the expression on the right of equation (11) will 
diverge to infinity with increasing t , either steadily or in the form of 
explosive oscillations, depending upon the sign and magnitude rela- 
tions of the members of 0,. In a certain particular case, however, 
namely when NX enters I, in suchia way that all the coefficients of the 
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terms in the 4;; « ©, are zero, the system will converge,.as before, to 
Si/ P i[E’, a, 8], if this lie within the region. This equilibrium point 
is highly special, since any infinitesimal divergence from the proper 
Initial conditions will lead to a non-zero coefficient of some term in a 
Ais € Q, and consequent explosion. We therefore term this equilibrium 
point an unstable one. We remark that the static analysis does not 
distinguish these essentially different types of equilibria. 

An especially interesting case is that in which either 0, is null 
or all the coefficients of terms in each j;; « 9, vanish, so that we do not 
have explosion, but ©; is non-null, so that some Aij = —1. If certain 
of these are multiple roots, and we find accordingly terms of the form 
aij;, %(—1)*, which diverge to infinity with t, we shall also suppose 
the coefficients of these zero, so that the z; all remain finite. In this 
ae we shall have asymptotically an expression for each z; of the 

orm 


2,—S./T;[E ,a, 8] + B(-1)', 


so that there are standing oscillations of constant amplitude in the 
steady state, whose amplitude is determined by the initial conditions. 
This will be a stable or an unstable equilibrium accordingly as we have 
had to assume zero values for the coefficients of particular terms in 
equation (11) to avoid explosion or not. 

We may now return briefly to the case where | — M| = 0, so 
that unity is a root of some of the T;. Here we may distinguish two 
contingencies: first, where the S; corresponding to every such T; van- 
ishes, and second, where this is not the case. The second case, as may 
easily be verified, is equivalent to the inconsistency of the equations 
(4) as discussed in the static analysis: here we shall have the par- 
ticular solution (9) for the z;; and (11) with this addend always di- 
verges to infinity, so that no equilibrium point of any sort can exist in 
T',—which accords with our previous conclusions. In the first of these 
cases, however, there is no particular solution to be added, and the 
appropriate discussion is exactly analogous to that for the case of 
Ai; = —1; if all coefficients of divergent terms be made zero, we shall 
obtain, asymptotically, z; = B,; + Bsi(—1)‘, where B,; and B,; depend 
on the initial conditions, except that B.; = 0 if @,; is null. We shall 
thus have, generally, oscillations of constant amplitude about a base- 
line determined by the initial conditions, which, in the case B,; = 0, 
leads to an arbitrary parameter in the expression for possible equi- 
librium points y in Ip; and since the number of T; with the root unity 
is equal to the nullity q of I — M, we find that there is a q-dimensional 
locus of such points in I’, , in consonance with the results of the static 
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considerations.” Which of these equilibria is*in fact attained is, of 
course, to bé'determined from the manner of entry of N into Ip. 
Let us defme a network-function (an N-function) in the follow- 
ing way: - 
(1). > (p; 2i)a? is an N-function, for any functions p; of peri- 
j=0 
od unity which’ do not vanish identically, and any constanta#1. 
(2). Any linear combination of functions of the form (1) is an 
N-function ; 7 


J Perl é 
(3). -N€x2)'+ > qjv? + K ( ) is an N-func- 


i PGF ley re) 
tion where N(x) ‘has the form (1) or (2), K is-a constant, the q; are 
uniperiodic functions which do not vanish identically, and 7 is zero or 
a positive integer. 

N-functions of the form (1) and (2) will be said to be of zero- 
order; those of the form (3) to be of order r. Moreover, we shall 
consider any:two N-functions to be equivalent if one arises out of the 
other by substituting any functions of period unity which do not van- 
ish identically for the p; and q;. The excitation in N within a given 
region I’, as:a. function of time, as given by equation (11), is an N- 
function and in normal form: we shall call it or any equivalent net- 
work function the characteristic function of I’) . Evidently, any char- 
acteristic function of I’, will serve equally well to specify the excita- 
tion-function of Nin I). 

Given these definitions, we are now in a position to state and 
prove a theorem which complements the foregoing results by a partial 
solution of the inverse problem, which is, to determine conditions un- 
der which a given set of excitation functions can be realized by a suit- 
able finite nerve-fiber network. We shall have, in fact, the 
THEOREM. 

Let a given P-space be partitioned into regions by planes perpen- 
dicular to the.axes in any desired way, and let a set of P functions be 
specified for each region, one for each coordinate axis: then, that 
there may exist a finite network N , with some pattern of applied ex- 
ternal stimulation, and having some set of P third-order synapses, 
such that the course of excitation at each such synapse when the sys- 
tem is in any of the regions I'p is given by the function specified for 
the corresponding coordinate axis in I’, , it is sufficient that the fol- 
lowing conditions be fulfilled: 


(A). Each of the specified excitation-functions must be an N- 
function. 
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(B). The given partitioning of P-space will define.a set of mul- 
tipliers analogous to those we have used above, though not 
necessarily univocally; and throughout every given region 
I, these multipliers will have a single value-distribution 
Tp1, Np2. Moreover, for every such distribution ap. , mp2 
such that no pair a; , B; are simultaneously zero, there is a 
corresponding region I). Then the condition is that, for 
some admissible set of multipliers, every region I’, whose 
specified excitation functions are not all constant must 
have at least one pair of multipliers a, , B, i Melcagel 


unity. 


In particular, given these conditions, it is possible to find a set of 
independent networks each of which consists of n simple circuits with 
one common synapse (we shall term these networks, which contain 
just one third-order synapse rosettes), such that N arises by run- 
ning chains from the centers of the rosettes to various designated 
points outside: but none back, so that the state of the whole network 
is determined by the states of the separate rosettes independently. We 
shall call networks of this kind canonical networks. 

In the proof, it will evidently be enough to construct a separate 
such NX for each dimension of P-space separately, since the N of the 
theorem is then the aggregate of these separate sub-networks. 

Now to every network function a 


T(x +1) 

N; (x) => Saya ary + Spo 2D P(e+1-—r)P(r+1) 
in normal form we may correlate a set of polynomials in E 
Yi" (E) = E(B — a4)" (EB — dig) ++» (EE — Qin) 


which differ among themselves only in a zero root of varying multi- 
plicity x. Let a suitable such polynomial be chosen, and denote the 
coefficient of H™ in it by O@m"(Ni), where n is the multiplicity of its 
zero root. 

Consider now a given region I’; , to which our hypothesis assigns 
a non-constant network-function N;, and let it have a set of charac- 
teristic multipliers of which one pair, say a; , 6; , corresponding to the 
limits A;, Y; are both unity. Now construct a rosette R in the follow- 
ing manner: if Y;‘(E) be of the n-th degree in E, R is to have 
nm — s circuits, Gn, Gna ,***»Gs, having respectively w,n—1,---,8 
fibers apiece. All the circuits are to have the same limits, namely A; 
and Y,— this can evidently be secured in every chain of sufficient 
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length, and here they are all longer than an arbitrary s—and the cir- 
cuit C, is to have an activity parameter A, = 6*n-»(Ni). We shall sup- 
pose that Be = oi + D mi;, 0; being the external stimulation at the cen- 


I 
ter of R; but o; and the w;; may be otherwise arbitrary. 
Now consider the course of activity in R. By the appropriate 
case of-equation (2) above, we find for this the equation 


Ee Wot (Ni) BP | y= W(E)y:=Be, (12) 
p= 


where y; is the excitation at the center of R. As before, we find the 
solution of this to be 


Vv. pis Tr 
Y¥:=> D dVijxvta,? + Dd prae* 
= 


j=1 k=>0 


P(E) T'(a +1) 
Derkeec| res ares (13) 


which becomes N;(¢t) when we set 
Bp,=6/([P(E)/(E—-1)7) aa. 


We shall suppose that rosettes R; have been constructed in this man- 
ner for every region J’; with a pair a; = $6; =1. 
Suppose now that 2;,, zi2 is the multiplier-distribution for I’; , 
and consider the function 
@,= vA (1 = a; ) (1 cae B;) IT a; B;. (14) 


4, JET G1 4,7ETi2 
Evidently, for the set of values for the multipliers 7;,; , rin , ; is unity, 


while for any other such distribution, it vanishes. ®; may be written 
out as a polynomial: 


P= > Wisx Oj1 Ajo *** Ajm Bra Bro +++ Ben» (15) 
d, 
in which every Wj;,= +1. 


The reader will easily see that if any two chains c, , c, have limits 
Ap, Yn; Aq, Yq, respectively, and corresponding multipliers a,, 8, and 
Oq, Pq, the result of putting them in series has the multipliers a, a, 
and #, 6,. It follows that taking each term on the left of equation 
(12), say Wax On One +++ Onm Bir Bro «++ Ben, We May construct a chain 
Cx. Whose lower multiplier is oj: on2 +++ Gam Bir -*+ Benj and we shall as- 
sign to it the activity parameter An = Wy», anda linn = 9. 

Now, taking the center of some one ®; of the constructed ro- 
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settes, say a synapse s;, and an arbitrary external synapse s;, con- 
nect a chain ¢, of this kind from s; to s; for each term on the right of 
equation (15). By the definition of @; , the excitation at s, will then 
be the same as at s; when the multipliers have the distribution qe 
%i; Otherwise it will vanish. Now, if we connect every rosette R; 
to s; in this manner, then whenever we are in some region I’; with a 
multiplier-distribution 2;;, ai2, we shall have yx = N;(t), which is 
the N-function required in the hypothesis. If, however, we have Yn = 
constant in J’; , we can have simply a single synapse s;, with a suit- 
able external stimulation, connected to s; in the same manner, instead 
of a rosette; and this will give the proper results at s,. Since we can 
make a construction of the above type for every dimension of P-space, 
the theorem follows. In particular, it will be noted that by this meth- 
od we may distribute equilibria of various types among the regions 
subject only to the conditions of the theorem. The necessity of the 
condition (A) of the theorem will be found evident. 


We may conclude by noting an immediate 


COROLLARY 

Given any finite network NX, it is possible to find a set of inde- 
pendent rosettes such that the excitation function of N for every re- 
gion is a linear combination of those of the rosettes—i.e., we can re- 
duce any network to a canonical network having the same excitation 
function. — 


In an intended sequel we shall consider the extension of results 
of the above type to networks governed by the two-factor excitation 
theories, instead of the present simplified linear model. We shall there 
develop the subject primarily from the standpoint of the inverse net- 
work problem, since it seems probable that it is here that the most 
fruitful and practically useful results are likely to be obtained. 

In conclusion I wish to express my appreciation to Dr. A. 8. 
Householder for his perspicacious counsel and criticisms. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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This first part of the paper deals with the theory of organic form. 
By introducing a formal principle, it is possible to develop a mathemati- 
cal theory of shapes of organisms, which covers both animals and plants. 


The development of Mathematical Biophysics has made consid- 
erable strides in the last years. Not only is the theoretical edifice wall 
erected and systematically developed, but a very large number of ex- 
perimental facts ranging from cell respiration to experimental aesthe- 
tics have been adequately quantitatively represented and in some cases 
even predicted by the theory. 

Naturally, as in every science, the more a theory is developed the 
more new problems appear in it. Also the more experimental facts 
have been explained, the more new experimental facts call for an ex- 
planation. The progress of any science may almost be measured by the 
width of the new horizon it uncovers to us. 

Whether there is such a thing as a possible exhaustion of a sci- 
ence, a stage in which everything is discovered and nothing remains 
to be done, is more a metaphysical question, and we shall not discuss 
it here. Premature worries in regard to the possibilities of explain- 
ing some of the yet unexplained phenomena, and anxieties lest such 
an explanation should turn out to be impossible, should not take place 
with scientists. Whenever such worries and anxieties have taken 
place, they usually were rather ridiculed by further developments by 
investigators who went ahead heedless of such worries. 

Nevertheless, at every stage of development of a science, when 
we have before us a vista of unsolved problems, it is useful, and per- 
haps even necessary to consider, so to say, the general plan of cam- 
paign in attacking those problems. We may not be expected to tell 
actually how to solve those problems, for that would practically 
amount to a solution. But we should be able to see at least in prin- 
ciple how those problems may be solved practically. 

If we survey from this point of view the field of mathematical 
biophysics, we find that while in some directions the method of further 
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conquest is already indicated by merely following and extending pres- 
ent techniques, in some other directions the adoption in the future of 
new methods and techniques is strongly indicated. 

The first case is exemplified mainly by the mathematical biophy- 
sics of the central nervous system. We may have to modify somewhat 
the fundamental equations. But in principle the rather formal method 
used hitherto, namely the study of geometrically more and more com- 
plex circuits, to which the formal, phenomenological equations are 
applied, still offers almost unlimited possibilities. 

' The second case is exemplified by the mathematical biophysics of 
the cell, especially of cell aggregates and of the organism as a whole. 
While by the introduction of drastic simplifications, a remarkable 
progress has been made in the dynamics of cell division, as well as in 
the theory of cell aggregates, yet an impartial survey of the situation 
shows that it is pretty difficult to even conceive in principle how the 
present methods are going to be applied practically to some more 
complex cases. We have outlined a couple of possible approaches to 
the problem of the form of-animals. In principle the development of 
those suggested approaches should well lead us to the solution of the 
problem. Practically, however, as it is easy to see even now, we shall 
soon run into almost insuperable difficulties if we try to extend even 
the present approximate treatment of the metabolic forces to such 
complex shapes as are offered by various organisms. Any oversimpli- 
fication of the problem must have a limit, beyond which the problem 
becomes completely distorted and unreal and a further simplification 
of the present “approximation method” would likely exceed that limit. 
Something else is to be done. 

Difficulties of a similar nature but perhaps even more severe in 
extent, are met when we remember.that hitherto mathematical bio- 
physics has studied only a very, very small section of the organic 
world. The interaction between the organisms and the surrounding 
environment has been taken into account from the very beginning. 
This interaction forms often a very essential feature of the theory, 
as for instance in the case of cell division. Yet the properties of the 
environment considered have been oversimplified almost to the limit. 
Only a very small fraction of the environment of an organism is in- 
organic. The largest part of that environment is formed by other or- 
ganisms. Almost every organism depends for its existence on the 
presence of other organisms, with the exception perhaps of some sim- 
plest unicellular organisms. With this exception the concept of an or- 
ganism or even of a number of organisms of the same kind, in the ab- 
sence of other kinds of organisms, is a fiction. 
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Again we may say, or at least we may hope, that in the final 
analysis all these complex interactions are reducible to physico-chem- 
ical metabolic forces of tremendous complexities. But that very com- 
plexity makes it again rather doubtful, whether practically we may 
ete anything by such a reduction. Again something else is to be 

one. 

Situations of that nature are not unfamiliar to the physicist. The 
latter does not doubt that even the most complex mechanical or elec- 
tromechanical phenomena are ultimately reducible to activities of in- 
dividual atoms. Nevertheless in studying for instance a complex elec- 
tric circuit, a physicist does not fall back on the equations of the elec- 
tron theory, but uses some general formal principles in his computa- 
tions. It is quite true that at present almost all such formal principles 
can be shown directly to be deducible from atomic theory. But there 
was a time when this was not known, and yet those formal principles 
were used just as safely. We may be reminded by way of example of 
Fourier’s theory of heat conduction, Ohm’s law, etc. 

We therefore propose to introduce in mathematical biology a few 
such formal principles, which may be used to advantage in the devel- 
opment of a theory of more complex biological phenomena, regardless 
at present of their reducibility to the principles used hitherto. It must 
be kept in mind, that the adoption of such formal principles does in 
no way compete with the use of the former ones, which are more 
physical and less formal in nature. The two sets are quite compatible. 
The introduction of the formal principles opens however new hori- 
zons, hitherto unseen. 

It may be worth noticing, that the mathematical biophysics of 
the central nervous system, which offers, as we remarked, almost un- 
limited possibilities for further developments with the present meth- 
od, is essentially based on a few formal phenomenological principles. 


I. Organic Form 

Whether a single cell or a complex metazoan or plant, every or- 
ganism is a metabolizing system. But whereas in a single cell meta- 
bolism is perhaps the most outstanding thing, complex metazoan or- 
ganisms are usually characterized by other essential physico-chemical 
attributes. One of the most important of those, which occur in single 
cells even, is motility. Another important attribute is the nervous co- 
ordination and integration. The velocity of locomotion of any or- 
ganism varies largely under different circumstances. We may how- 
ever consider every mobile metazoan as endowed with an average 
velocity of locomotion, v , characteristic of the species. Similarly the 
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amount of nervous coordination and integration varies largely from 
case to case. But again we may consider average values, charac- 
teristic of the species. This average amount of coordination will be 
some increasing function of the amount of nerve tissue per unit 
mass of the organism. 

With the exception of some protozoa, the locomotion of all or- 
ganisms is produced by a system of levers. In such cases as snakes or 
fishes the body as a whole is used as a lever. In other cases the ex- 
tremities are such systems of levers. 

Although quite a considerable amount of research has been done 
on the kinematics of movements of the vertebrates and of some other 
organisms, (Bethe 1930), a systematic theory of what may be called 
lever-propelled metabolizing systems is yet completely lacking. An 
elaboration of such a systematic abstract theory, as a preliminary to 
actual biological applications, is strongly indicated. 

Given a metabolizing system of mass M, average metabolism q 
ergs gm™ per unit mass, average amount of nervous tissue nM gm., 
the average velocity of locomotion v will be determined by the partic- 
ular structure of the systems of levers which constitute the locomotory 
system. If we have a systematic classification of such systems of 
levers and of their properties, we can calculate for a given lever sys- 
tem, the value of v in terms of M , q and 7. In general we may expect 
that for given M, q and n the same value of v may be obtained by 
means of several lever structures. 


In developing a systematic theory of lever propelled metaboliz- 
ing systems, we shall amongst other things classify those systems by 
their degree of complexity. Provisionally we may consider the com- 
plexity as measured by the number of single rigid levers of the sys- 
tem. A more detailed theory should, however, take into consideration 
also the different types of joints between the levers. A joint with two 
degrees of freedom will be considered as more complex than a joint 
with one degree of freedom. In such a case the complexity may be 
measured by the number of kinematic degrees of freedom of the whole 
system. 

We shall now introduce the following formal principle. 

A lever propelled metabolizing system, characterized by given 
values of M,q,nand v has the simplest lever structure which for the 
given values of M , q and n gives the required value of v . 

Since the lever structure of an organism determines to a large . 
extent its shape, we thus may describe the shape of an organism in 
terms of M,q,nandv. 

Not much can be done with the above principle until the actual 
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theory and systematics of lever propelled systems is sufficiently devel- 
oped. In order to show however in what direction such studies should 
te place, we shall give here a few examples, meant as illutyations 
only. 

Let us first consider the simplest case when the system is of a 
cylindroid shape, flexible along the axis and can itself be used as a 
folding lever system. In this case which in nature is exemplified by 
snakes, the locomotion proceeds essentially by alternate folding and 
unfolding. During the folding the “front” end (Figure la) remains 


—> —> 
a b c 
FIGURE 1 


in fixed contact with the ground, during the unfolding the “rear” end 
does it. (Figure 1b). In the locomotion of a snake actually this fold- 
ing is a double folding (Figure 1c), but this does not change thing 
in principle. : 

This type of movement requires relatively little coordination and 
therefore a relatively small n will be sufficient. We may set as require- 
ment n > 7,, where n, is the minimum n required for such types of 
movements. 

For simplicity (since this is an illustration only) let the angle of 
folding be constant. Then, denoting by » the frequency of folding and 
unfolding, by L the total length of the system, and by a a coefficient, 
we have 

yv=aL. (1) 


The energy necessary to maintain the system in continuous motion 
with the average velocity v must be supplied by the total metabolism 
qM . That energy is essentially composed of two parts. The first part 
is the overcoming of the frictional resistance of the ground. In dry 
friction the resistance is independent on the velocity. In that case the 
work done per unit time is proportional to v. If the resistance is 
proportional to v , then the work done per unit time to overcome that 
resistance is proportional to v?. Moreover, for dry resistance it will 
be proportional to M , for the friction will increace with M. The sec- 
ond part of the energy is that required to actually produce contrac- 
tions of the muscles involved, and is per unit time proportional to the 
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mass of the muscles and to vy. Since in the case considered the whole 
body participates in the locomotion the mass of the muscles may be 
put equal to M. 

Denoting by c and f two more coefficients and considering that 
practically the whole metabolism contributes to the locomotion, we 
have for the case of dry resistance 


qM =cMv + frM , (2) 
or 
eat (3) 
Equations (1) and (3) give 
Aw ae (4) 


The length L for a given shape varies as M‘?. Putting for sim- 
plicity the average density 6 of the system as equal to 1 (which it 
practically always is) we may write 


L=aM’, (5) 


a being a form factor. Considering the shape of the organism as a 
very oblong ellipsoid of revolution with major semiaxis 7, and minor 
T2, and putting 


UY owe 
me (6) 
we have 
4 
3 ur, To" == M 
or because of equation (6) 
4 uri 
3% prea (7) 
Hence 
pave eee One 
L=2r,=2 [mn (8) 
Equations (5) and (8) give 
_ */ 3A? ois OE 
re es i=. (9) 


From equations (4) and (5) we have: 
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__ a(q— cv) «o M2 
f i? 


v (10) 


or 
fv 


ne A118 e=ren): (11) 


Hence, if for such a system, M , q and v are given with n > 1, 
then a is completely determined. But a determines by means of 
equation (9) and 4 determines L by means of equation (8). Hence 
the gross form of the organism is determined. 

The quantity a, and hence also 4 increase with v. For 

vas (12) 
a becomes infinite. Hence the value (12) of v is the greatest that can 
be obtained for any M and q, with the particular form of the organ- 
ism. 

If 

ee (13) 
C 
then the organism must represent a different system of levers, and 
must possess lever attachments in the form of extremities. We shall 
presently see, that introduction of extremities does actually raise the 
limit of v over that given by equation (12). 

In principle we may consider any number of pairs of extremities. 
For simplicity consider them all of the same length 1. They may work 
either alternately or synchronously and the locomotion again consists 
in a folding and unfolding of each pair, the “rear” member of the 
pair pushing against the ground in the unfolding process, the “front” 
member during the folding. 

Let us first consider the purely theoretical case that the extrem- 
ities are rigid levers of length /, that is without knee-joints, etc. Then 
denoting by a, a constant, we have, similarly to equation (1): 


v=a yl, (14) 


regardless of the number of pairs of extremities. Hence, it would at 
first appear as if our fundamental formal principle would allow only 
organisms with only one pair of extremities. This however is not so, if 
we consider the role of the quantity n. Both locomotion and quiet stand- 
ing are unstable processes for only one pair of extremities. The main- 
tenance of standing posture with only two extremities requires a very 
complex mechanism of nervous coordination and hence a sufficiently 
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large n. Let us denote by 7, the smallest value of n which will pro- 
vide for the necessary coordination. As a measure of n, we may per- 
haps take the total number of muscular movements, which is involved 
on the average per unit time in order to maintain an erect posture on 
two extremities. 

If n = nm then the organism may have only two extremities. This 
does not mean that it necessarily will have. For if inequality (13) is 
not satisfied, the minimum number may still be zero. But if (13) is 
satisfied, and n < n,, then the organism must have at least two pairs 
of extremities of the rigid type assumed above. 


A standing posture on four extremities requires much less co- 
ordination, although mechanically it is not quite stable, considering 
that the extremities are not rigidly joined to the body. A certain 
minimum amount of muscular coordination is still necessary to main- 
tain a standing position. Let to that minimum correspond a value % . 
Then we find four extremities only within the range n, > n > n,. For 
Nn, > Nn > nN, we find six extremities, etc. 

There is another factor which plays a role in the determination 
of the number of extremities. That is the size of the organism relative 
to the average “roughness” of the ground on which it is moving 
around. Even a human uses all four extremities when climbing over 
a rough mountainous terrain, when the average size of the “bumps” 
becomes of the same order of magnitude as that of our bodies. For 
very rough terrain we have to use our fingers, which are additional 
levers, to take hold of the unevennesses. Hence a system of even four 
rigid extremities will not be sufficient in such a case, and it is neces- 
sary to introduce additional levers either by adding more joints at the 
ends of the extremities, or by adding more extremities. 

From kinematic considerations it may be readily seen, that in 
the first instance at least two “fingers” have to be added to each ex- 
tremity, to provide for sufficient hold. Hence, the second case intro- 
duces less additional joints. The smaller number of joints the smaller 
may be n. Hence for small values of 1 the stabilization will take 
place by addition of new extremities. And since the relative roughness 
of the ground in general increases as the size of the organism de- 
creases, we shall find 6 or 8 extremities for small values of both M 
and m. Such is the case with insects. When M and 7 are large, only 
four extremities may be sufficient but they will be provided with “sec- 
ond order” extremities. 


Now we may consider the actual expression for v in terms of M 
_and q. In addition to the overcoming of resistance to the motion and 
to supplying the energy of muscular contraction, the organism mov- 
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ing by means of rigid extremities must spend with every step some 
energy on lifting itself (Figure 2). This energy is not recovered in 


—> 


FIGURE 2 


the downward motion, or at least not completely recovered, due to 
friction etc. Per unit time it is equal to 


br» lM (15) 


b being a coefficient. 

For the work to overcome the resistance of the ground we shall 
have the following expression; we may consider that the resistance 
is exercised everytime the extremity, during the forward motion 
comes in contact with the ground. That resistance force is roughly 
proportional to the velocity v, and hence the work of resistance is 
proportional to v?. But this work is done only intermittently, the fre- 
quency of this intermittent work being ». Hence the total work is 
roughly proportional to 


M yv?. (16) 


Denoting by M, the mass of all the extremities and remembering that 
now only those are involved in the locomotion, (which is of course 
only a crude approximation) we have now instead of equation (2): 


qM.=bryIM+¢,Mrv+f.»M., (17) 


where c, and f,; are different from c and f. 
Solving equation (17) for » and introducing it into equation (14) 
we find 
Bef sit Sur Oar (18) 
bM1 ela c,Mv? a5 f.M. 


Introducing a “form factor” § for the extremities, we may put 


r= pio, (19) 


Vv 


and equation (18) now gives, with 6’ = 1/f*: 


dnBiqi (20) 


°=FMI + ¢,Mv? + BRE’ 
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which may be written thus: 
a,p’qlt — B'fyvlP — bMvl = ¢,Mv. (21) 


Equation (21) determines / for a given set of M,q,» and f’. It may 
be shown that equation (21) has always one and only one positive 
root, /,, which represents our solution, the negative roots being dis- 
carded as physically meaningless. The form-factor f’ is determined by 
mechanical considerations. For a given M, the length of the extrem- 
ities is limited by their mechanical strength needed to support the to- 
tal mass M of the organism. As v increases, J, increases. Hence for 
any M , q and v we may determine a corresponding / = 1,. Therefore 
if inequality (13) is satisfied, crawling is impossible, but running may 
be possible. The only limitation imposed now is the circumstance, that 
for a given M , 1 has an upper limit /,(M) for purely geometric rea- 
sons, since an organism cannot be “all legs”. If inequality (13) holds, 
and if for that v the necessary / is greater than I, , then neither crawl- 
ing nor running is possible. : 

Consider however that each extremity has a “knee joint” and 
can be folded. In that case the term with b in equation (17) and the 
following ones may be reduced or eliminated, for there is no need of 
lifting the organism (Figure 3) as much as before. The result will 


—_> 


FIGURE 3 


be a decrease of J for the same values of M, gq and v. Now we may 
have l, < 1, for v > q/c, and while crawling is impossible, running 
is possible. 

If the constants in our equations, which have to be determined 
by a more exact elaboration of the theory of lever propelled systems, 
are such, that J, > 1, whenever b > 0, then for such values of vy which 
satisfy inequality (13) any locomotion will be possible only if the ex- 
tremities have a knee-joint. That this is so, will be demonstrated in 
a subsequent paper. 

Making more than one knee-joint will not improve anything, and 
therefore, according to our fundamental principle, such a structure 
does not occur. 

In a similar way we may treat other types of locomotion on the 
ground. For instance a system based on a structure without extrem- 
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ities may propel itself by successive folding and refolding in the verti- 
cal plane, as in the case of some caterpillars. Here the term in v? in 
equation (2) vanishes, but a term similar to that given by (15) ap- 
pears. A limitation is also imposed by considerations of mass and 
rigidity. It will be readily found, that for a sufficiently large mass 
the system will not be rigid enough. The necessary inequalities may 
be readily derived. But even without deriving them it may be seen 
that such type of locomotion will occur only in small animals. 

All the above is given only as an illustration, pending the neces- 
sary development of a systematic theory of lever propelled systems. 
But our illustration shows, how by combining our fundamental prin- 
ciple with purely mechanical and physical considerations, we may de- 
termine the given form of an organism in a unique way for a given 
set of M,q,andv,n. 

Locomotion in water and air is treated essentially in a similar 
way. Consider for instance locomotion in water by the use of the 
whole body as a lever (some fishes). We have again, denoting the 
body length by L and by a, a constant 


Y= O7L . (22) 


The resistance now varies approximately as the linear dimensions 
or M2, Therefore we have now 


qM = forM ae C,M13y2 . (23) 
This gives 
_ QML — c.M*Lv? 


4 
7M Sea 


and equation (24) determines L in terms of M,qandv. 

The treatment of flying is somewhat more complicated, requiring 
aerodynamical considerations. But in principle it proceeds in a simi- 
lar way. 

It may be readily seen that less energy is required to keep a small 
organism afloat in the air than a big one. If an organism is falling 
through air, the air resistance for not too large velocities is approxi- 
mately proportional to rv , where r is the average linear dimension. 
(Stoke’s law). The force F of weight, pulling it down varies as oe 
Hence v ~ 72. The work done per second by the force of gravity is 
proportional to Fv ~ r°, while the total metabolic energy available is 
proportional to 7°q. Hence, from 7° = 7°q we find that 


Gor, (25) 


It may be therefore expected that the aerodynamic efficiency of any 
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mechanism used for flying, may be much smaller for small organisms, 
than for large ones. Hence also the degree of nervous coordination 
and therefore n, will be much larger for larger organisms. 

We may next consider organisms, that can propel themselves in 
two or three different. media, e.g. ground, water and air, with corre- 
sponding average velocities v,, Vw», V2 and determine in the above 
manner-the sizes and number of their legs and wings. As we have 
seen for n > n., only two legs will be found. With a large M, as we 
just have seen, a large 1 is required for flying. If that value of n = n,; 
is always larger than n,, then all large flying animals will have only 
one pair of extremities for running (birds). For small M , flying does 
not require such a large n , and the number of extremities for ground 
locomotion will be larger, as determined by previous considerations. 

Thus far we have considered only animals. With some excep- 
tions, plants may be considered as organisms with v = n= 0. The 
problem of form may however be treated also in this case on the ba- 
sis of similar principles. 

First of all, the absence of locomotion and therefore the impos- 
sibility of running after the food, will require a much larger specific 
surface of the organism through which food can be absorbed from the 
surrounding medium. Hence the general “branched” character of 
plants. As a further illustration, we shall outline a possible approach 
to the form of plants. 


Let 1, be the length of the main trunk, 7, its radius. Let 1 and r 
be the average length and radii of the branches. Let the total number 
of branches be n , and the average density of the plant 6. Then 


M = 3(lor.? + ult"). . (26) 


Since gM is proportional to the total number of leaves, it is therefore 
approximately proportional to the area of the branches. Hence 

qM = kuir , ; (27) 
k being a coefficient. | 


From mechanical considerations, the length of a branch cannot 
exceed a certain value, determined by its radius, lest the branch break 


L=f(r,6). (28) 


Since the smaller 7 , and the greater 1, the greater for a given M the 
specific surface of the tree, if the latter has a maximum, the relation 
(28) must hold. 


Similarly 
. l=fo(%,6,M). (29) 
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Moreover, the total flow of metabolites through the trunk is limited 
by r), and the diffusion properties of the material. Considering that 
7. is always just such as to allow the necessary flow, we have for the 
trunk 


1) =f: (qM , 0), (30) 
and for each branch va 
\ 
r=fi(+—, 0), 5 (31) 


Altogether we have six equations, to determine the six quantities J, , 
To,n,l,rand 6 in terms of M and q. But the first five quantities al- 
ready largely determine the gross shape of the plant. Small J, and 7, 
and large n, give us a bush. Small J,, large 7,, small #-and large 1 
give us the form of Figure 4a. Large 1,,7.,n, land 7., give:the shape 


eile 


—. 


a b 
FIGURE 4 


of Figure 4b, while the large 1, , 7, and ” , but small / and r give some- 
thing like Figure 4c. By using plausible expressions for f(r, 4), fo 
and f, those different forms are actually obtained as functions of M 
and q. This will be discussed in detail in a subsequent paper. 

Instead of equal branches we may introduce a distribution func- 
tion n(1) , such that f2n(l)dl=n, and determine n(l) from some ex- 
tremum requirement. iat es 

The shapes of leaves may possibly be related to their.mass and 
metabolic rate. Thus consider a leaf as shown on Figure 5a. Let its 
mass and hence its size, increase. The distance between the.veins be- 
comes greater, until it eventually becomes too great for the metabo- 
lites to reach the remotest points. Hence 5a becomes impossible, 
while shape 5b will be found. The same effect is obtained by increas- 
ing q for a constant mass. ho 


46 GENERAL BIOLOGY 


FIGURE 5 


Different shades of green of a leaf are determined by the chloro- 
phyll content, and are therefore also linked to g. Combining all these 
considerations together, we see a way of establishing a relation be- 
tween M and q, the general shape of a plant, and the size, shape and 
shade of its leaves, and through 6, the mechanical properties of the 
wood. 

The gist of the whole consideration is that we express the essen- 
tially discontinuous properties of the organic form (no extremities, 
two extremities, six extremities) by a set of continuously varying 
parameters. We do not need to limit ourselves to M,q,v and n, and 
we may apply the same procedure to cases, when there are a larger 
number of determining parameters. Thus in addition to the average 
velocity of locomotion v we may consider the maximum acceleration v 
of which an animal is capable and which play an important role in 
its life. This may lead to a theory of jumping, involving considera- 
tions of unequal hind and front extremities. 

The structure of different sense organs, with which an organism 
is endowed, may also be treated in a similar way. We shall consider 
as continuously varying determining parameters such things as acuity 
of a given sense, range of perception, etc. and by developing first an 
abstract theory of perceiving mechanisms, infer from the values of 
the determining parameters to the geometrical and physical structure 
of the sense organ. 


As an important parameter we may introduce the degree of dif- 
ferentiation of an organism. Tentatively this may be measured by 
the average maximum amount M, of M that can be removed without 
destroying the organism as such. We may for instance put for the 
degree of differentiation D 
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M 
D= Me (32) 
While all the above is at this stage meant only as an illustration, 
even in this crude form the theory suggests a number of new ex- 
perimental approaches. The relations between M, q, v, » and the 
length of extremities can be in principle verified experimentally. The 
same holds of equation (4) and (24). The general conclusions even 
of this sketchy theory show some overall agreement with actual facts. 
Thus we see the reason why large flying animals (birds) use only two 
extremities for running, while smaller ones have six or more. The 
appearance of four extremities for running in higher developed ani- 
mals and its reduction to two in the highest, becomes also intelligible. 
A relation between length of extremities and size of the nervous sys- 
tem suggested twice before from entirely different considerations, 
now appears in new light. Much is still to be travelled, but the road 
seems to be clearly indicated. 
The author is indebted to Dr. Alston S. Householder and to Dr. 
Herbert D. Landahl for reading the manuscript and valuable com- 
ments. 
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